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Abstract. Using the £ functional equation and the Hadamard product, an 
analytical expression for the sum of the reciprocal of the £ zeros is established. 
We then demonstrate that on the critical line, [£| is convex, and that in the 
region < 5R(s) < 0.5, |£| has a negative slope, given the RH. In each case, 
analytical formulae are established, and numerical examples are presented to 
validate these formulae. 



1. Introduction 
The Zeta function was first introduced by Euler [4] and is defined by: 



(1.1) C(s) = 



oo 1 

E- 



Riemann [6] extended this function to the complex plane, for 3?(s) > 0, mero- 
morphic on all of C, and analytic except at the point s = 1 which corresponds to 
a simple pole: 

^ ^(^r^yE — 

The Zcta function satisfies the functional equation: 

(1.3) C0O = x(*)C(i-«) 

where: 

(1.4) x(s) = 2 s t: s - 1 sin(.57rs)r(l - s) 

Let s* be a zero root of £(s) such that |C( S *)I = 0. There are a number of "trivial 
zeros" at s = —2,-4.... All the other (non-trivial) zeros must lie in the critical 
strip < 5R(s) < 1. Clearly, equation |1.3| is satisfied if 5R(s*) = .5. This is the 
Riemann Hypothesis - the nontrivial zeros of £(s) have a real part equal to 0.5. 
The RH is proved if one can demonstrate that there are no other zero roots - that 
is, all the roots fall on the critical line. 

The approach taken in this paper is to first derive analytical expressions for the 
derivatives of \x\', this occurs in Section 3. Then in Section 4, we derive analytical 
expressions for the derivatives of |£| by exploiting the functional equation and the 
logarithmic derivative of the Euler product using the Hadamard product. An an- 
alytic expression for the sum of the reciprocal of the C zeros is established, which 
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permits a statistically based conjecture about the truth of the RH. We show that on 
the critical line, |£| is convex and that in the region < 5R(s) < 0.5 and t > 2tt + A, 
\(\ has a negative slope, given the RH. At each stage we demonstrate the validity 
of the expressions derived using numerical examples. 

2. Analysis 

We follow Riemann's notation: s = a + it. In this work, we take 3(s) = t as 
fixed, following Saidak [7,8]. 

Lemma 2.1. The non trivial zeros of the C function lie on the critical line, or 
occur in pairs equidistant from the ciritical line. 



This follows directly from the functional equation 1.3 - if C(s) is zero, then so is 
C(l — s). Thus the complex zeros of £ cither lie on the critical line, or are symmetric 
about it in the strip < 5ft(s) < 1 Since such zeros occur in pairs, any violation of 
the RH can be investigated by considering just one side of the critical line. In this 
analysis, we consider the area !ft(s) < 0.5. 

The following properties of the £ function are used throughout. On the critical 
line s = .5 + it, 

nc(s) = 5RC(1 - s) 

9fC(«) = -9C(l-«) 
argC(s) = -argC(l-s) 

In addition, since we are dealing with absolute values |/(s)|, for s = a + it, we 
note: 

ln/(s) =ln|/(s)| + iarg/(s) 

and hence: 

(2D \_mi_ n (m 

[ } \fis)\ *\f(a) 

where f'(s) implies d ^ evaluated at s. 

We will be dealing with a particular region of the critical strip, which we define 
as s e S : < SR(a) < 1 and 9f(s) > 2n + A. 

3. Analysis - \\\ function 
The Zeta function satisfies the functional equation: 

(3.1) C00=x(s)C(l-s) 
and hence: 

(3.2) IC(«)I = |X(*)I|C(1-*)| 
where: 

(3.3) x(s) = 2 s tt s " 1 sin(.57rs)r(l - s) 

In this section, we derive the properties of |x( s )| which are required in ascertain- 
ing the properties of |C( S )I- I n particular, we consider the analysis of |x(s)| t=to - 
the slice of |x( s )l when t is fixed at to. 

A graph of |x(s)| , s € S is shown in Figure 1 - the second contour depicts the 
critical line 3?(s) = .5. 

The following properties of |x( s )| are established: 



TOWARDS A STATISTICAL PROOF OF THE RIEMANN HYPOTHESIS 




Figure 1. 3D plot of \ x (s)\, R(s) = 0...1,3(s) = 2tt..200. 
Lemma 3.1. For s G S, |x( s )l * s continuous. 



From equation 3.3 for s G S, x(s) is a holomorphic function, and thus |x( s ) 
continuous. 

d 



Lemma 3.2. For sGS,£ |x(s)U = |x(*)| M**) - - a))) 

For any i = t : 

'| = 2 CT ln(2) 



_9_ 



i 



da 



|sin(.57rs) | 



■k"' 1 ln(7r) 

.57T sin(.57r cr) cos(.57rcr) (cosh 2 (.57rfo) — sinh 2 (.57rf )) 



|sin(.57rs)| 



9(7 



|r(i - S ) 



.257r sin(7r tr) 
|sin(.57rs)| 

-|r(i- s )| -*)) 



where is the digamma function, and the last equation follows from [2~T 
Taking the modulus of equation 3.3 we have: 

(3.4) |x(s)| = |x(<7 + **o)| = |2 S | k 8 - 1 ! |sin(.5™)| |r(l - s)| 
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Hence the derivative of \x{ s )\ with respect to cr, while keeping t = to, is: 

d . / \ i i / \ i /i / n \ .257rsin(7rcr) „,,-,,., ss , 

ozlx(*)lt=t = IxWl ( M2tt) + , \ a - ft(*(l - a)) ) 
w |sin(.57rs)| 

(3.5) = \ X (s)\Hs) 

The first term in (j)(s), ln(27r), equals 1.837877067. The second term is positive 
in the critical strip, but quickly becomes insignificant - at t = 2ir this term is less 
than 0.0001, and by the first root (s = .5 + il4.134725) it is less than 4 lO" 10 . The 
last term 3^(^(1 — s)) is the real part of the digamma function (which is asymptotic 
with the natural logarithm function) and is positive in the critical strip for t > 1.1. 
For A = 0.02, 4>{s) < for t > 2tt + A. Thus, for s G S, and t > 2tt + A: 

(3.6) ~ \x(s)\ f=tQ = \ X (s)\ (M2tt) ~ - *))) 



Equation |3.6 
between equation 



da 

works extremely well. Using Maple, a numerical comparison 



3.6 and the numerical derivative of |x( s )l evaluated over the 



critical strip at the first root showed a percentage deviation of less than 1.4 10 7 , 
with even lower deviations for a =/= .5. 

Lemma 3.3. For s G S , the minimum and maximum of \xi s )\t=t occurs at a = 1 
and cr = respectively. 



For s € S, x(s) in equation 3.5 is a holomorphic function with no zeros. Hence, 
the extreme points of |x( s )lt=t wm " ne on the boundary. 

Lemma 3.4. For s G S, \x(s)\ > 0. 

For s — a + ito, let: 

$(s) = J (sin 2 (.57rcr) cosh 2 (.57rf ) + cos 2 (.57rcr) sinh 2 (.57rt )) 

and thus: 

|x(a + i < )| = 2 CT ^- 1 |r(i- s )|<i>( s ) 

and using standard results [1]: 



\n-ih))\ = 



|r(i-ito))| - 

Hence for s G S, and from Lemma |3.3[ 



7rf 



sinh(7rfo) 

/ ; * 

to sinh(7r<o) 



lim \ x {l + it )\ < lim \x(e + it )\ < lim \x(Q + it )\ 

Zq — S-OO Zq — *oo Iq — s-oo 

(3.7) J^- < lim \x(v + *t )\<\[]f 

Thus, |x(s)| > 0. These limits apply even for relatively low values of to. Using 
Maple, a numerical comparison between ^t /27T and |x(0 + ito)| evaluated at the 
first root showed a percentage deviation of less than 0.7 10~ 7 . 

Lemma 3.5. For s G S, t > 2tt + A, |x( s )l t=to * s negative. 
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From Lemma [3. 2| 
d 



From Lemma 



3.4 



()u \x(s)\ t=t0 = \x(s)\ (ln(27r) - - a))) 

\x(s)\ > 0. Thus, for small values of t (t < 2n), \x(s)\ is 



da 

positive, dominated by the first term, ln(27r). However, for < er < 1, t > 2n + A, 
the last term 3?(W(1 — s)) dominates, and hence g^|x(s)| t=to < 0. This also 
establishes monotonicity. 

Lemma 3.6. If3t(s) = .5, then \ X (s)\ = 1. 



On the critical line, \((s)\ = |£(1 — s)|. From the functional equation 1.3 

IC(*)I = \x(s)\ |C(i-»)l 

(3-8) = |x( S )| |C(«)| 

and hence |x( s )l = 1- 

This result also follows from the definition of x(s). For any t = to and for 
o = 0.5 '. 

|2 S | = 2 a = V2 
Itt 5 - 1 ! = tt* 7 " 1 = I/Vtt 

|sin(.57rs)| = y sin 2 (.57T(r) cosh 2 (.57rio) + cos 2 (.57rcr) sinh 2 (.57rfo) 



/sin 2 (7r/4) cosh 2 (.57rf ) + cos 2 (tt/4) sinh 2 (.57rf ) 

|r(i-s)| = |r(.5-»to)| = VVcosh(^o) 

Thus from equation |3.3| 



(2/V)(.5cosh 2 (.57rf ) + .5 sinh 2 (.57rf ))7T 
cosh(7rio) 



|x(.5 + **o)l = \ 
= 1 

Hence, if 3?(s) — .5, then |x( s )| = 1- 
Lemma 3.7. ForseS,t> 2vr + A, |x0)|" = p^yp 



Differentiating equation 3.5 provides the second derivative of \x\' 



For even moderate levels of t, d 3 ^ 1 s ) _ q jjcncc 

(\x(s)\'f 



\x(s)\" = 1x001 ( dm ^ s) + H27T) - m(i a )fj 

derate levels of t, g^j-llrf) 

7 e'er 

(3.9) IxMl"- Ws)| 

Using Maple, a numerical comparison between equation |3.9| and the numerical 
derivative of equation |3.5| evaluated at o — .5,t = 200, showed a percentage 
deviation of less than 1.7 10 -7 , and even lower off the critical line. 

At ffi(s) = .5, |x(s)| = 1, and thus 

\x(s)\" = (\x(s)\') 2 
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3.9 



Lemma 3.8. For s G S, t> 2vr + A, \x(s)\ (k) = f^jl-i 

The fcth derivative of |x| can be derived by repeated differentiation of equation 
We note that the sign of the fcth derivative is (— l) k . Also, for 3?(s) = .5, 

i\0)l w = (lxO)l') fe 

4. Analysis - |(| function 

In this section, we derive the properties of |C( S )| 
Lemma 4.1. If s* is a zero root o/£(s), then |C( S *)I = 0- 

By construction. 

Lemma 4.2. For s£ S, 9ft(s) < 0.5, |C0)| > (CfT^s)! 

From Lem ma |3 .5[ for s £ S, \x( s )\t=t Q ^ s nionotonic with ^ |xO)lt=t < 0. 
From Lemma pT6T]y(0.5 + it )\ = 1. Thus for 3?(s) < 0.5, \x(s)\ > l.°Using 



equation 3.3 for 5R(s) < 0.5, and noting that |C(s)l = |C( S )I ; 



(4-1) |C0)l = 1x0)1 |c(i-«)| = 1x0)1 \C(i-s)\ 

> |C(T^)| 

Remark 4.3. If the strict equality held, we would have RH; however the equality 
holds at the zeros. 

Lemma 4.4. For ffi(s) < 0.5, |C(s)|' < if - s)\' > 0. 

At any point s to the left of the critical line, |C0)l can b e analytically evalu- 
ated as a function of the slope at the corresponding point 1 — s to the right of the 
critical line: 

CO) = xO)C(i-*) 
|C0)| = |xO)IIC(i-*)l 

(4.2) icwr = -ix^nca-^r+Ka-^iixwr 



|x(s)| and |£(1 — s)| are positive by construction. By Lemma 3.5 |x0)| is neg- 
ative always in the critical strip. Thus if — s)\' is positive, then |C0)| nas to 
be negative. 

Remark 4.5. Unfortunately, this does not help for the situation of interest - that is 
where - s)\' < 0. 

Lemma 4.6. For ft(s) = 0.5, \((s)\' < 

From equation |4.2[ 

(4.3) |C0)l = - 1x0)1 IC(i-*)l' + IC(i-*)l 1x0)1' 

At s = .5 + it, |C(s)| = |C(l-s)|, |C(l-s)|' = |C(s)|', and by Lemma 
1x0)1 = !■ Hence, 

|C0)l' = .5|C0)||x0)l' 

< 



3.G 
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since by Lemma 
(4.4) 



3.5 



|x(s)f < Thus at K(s) = 0.5, 
\C(s)\' 



\C(s)\ 



Using Maple, 9? «'(s)/C(s)), evaluated at s = .5 + 200i, gave -1.7302196, while 
■5|x(s)f gave -1.7302197. 



Lemma 4.7. 



Ep*(?) 



log(2) — .51og(7r) + 1 + .57 where p are the non trivial 
zeros of the Q function, and 7 is Euler's constant. 

For any s in C, the £ function is defined as: 

£(s) = .5s(s- l)7T' 5s r(.5s)C(s) 
The Hadamard product representation of £ is: 

p 

where p runs over all of the non-trivial zeros of C( s )> an d A and B are absolute 
constants. Taking the logarithmic derivative of £(s) we obtain: 



B 



E 



Rewriting in terms of £, we obtain: 



(4.5) 



= B — 



1 



1 



.51og(7r) — .5 



r 7 (.5g + i) 
r(.5s + i) 



E 



Evaluating equation 4.5 at s = gives: 

B = log(2) + .51og(7r) - 1 - .57 
where 7 = -r'(l) is Euler's constant. Thus B = -.0230957 
Since 



icmi J - x \ a*) )> 



4.4 



and 



4.5 



we can compare equations 

(s) = 0.5. The term —3=- tends to zero as t gets large, and using Lemma 
note that: 

\x[_ 
\x\ 



evaluated at 
we 



3.2 



(4.6) 



ln(7r) - SR(^(.5s + 1)) = ln(27r) - 3?(*(1 - s)) 



Thus, if these two equations are to have the same value at 3?(s) = 0.5, then: 



(4.7) 



p 



1 



= 3? 



= —B 



Evaluating the sum over the first 100000 non trivial £ zeros for s = .5+it for a range 
of t, t > 200, generated a constant value of 0.023073645, which is approximately 
-B0 

Lemma 4.8. For R(s) = 0.5, |C(s)|" > 0. 

^The difference comes about because of the summation over only 100000 terms, instead of 
infinity. We noted also that the sum was over both positive and negative values of the zeros - and 
indeed in Hadamard's work the three regions were defined in terms of the absolute value of p. In 
addition, the analysis is not valid for a situation where s is equal to ft, a ( zero, since |£(s)| is 
not analytic at that point. However it is valid in the limit. 
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Differentiating equation |4.5 



(4.8) 



C'(s) fC(s) 



as) 

Note that lim t _, 
»(«) = 0.5, 5R (E P 

(4.9) 

and equivalcntly: 
(4.10) 



C(«) 



1 



1 



*'(. 5s + l)-£ 



1 



(s - p) : 



00 (s-1) 
1 



(S-P) 2 



(a -1)3 

= and limt^oo ^'(s) = 0. In addition, evaluated at 
= SR(Jf (2, s)) < 0. Thus we have: 



if(2,s) > 



Of) _ 



C(«) 

_ f\C(s)\' 

\C(s)\ \\t(s)\ 



-K(K(2,8)) 



Again, we check to see that this is correct numerically. Using s 



200z, 



CM 



CM 



1.4511544. Evaluating £) p f ^j^p J over 100000 zeros gen- 
erated K(2, s) = -1.4511540. 

Remark 4.9. The evaluation of K(2,s) = J2 P ^ ( (s-p) 2 ) ' ^ s on ly a function of t, 
since 5R(s) = 5i(p). Defining K(k, s) — Sft J , it is clear that for k odd and 

R(s) = .5, 5R(iT(fe,s)) =0. 

Lemma 4.10. For < 3?(s) < 0.5 and f > 27r + A, and given the RH, |C(s)|' < 0. 
Applying Lemma |4.7| and equation |4.6| to equation |4.5| gives : 



IC(*)I 



where as before we ignore the term SR^y which becomes negligible for large t. 
Rewriting as: 

'lx(*)l' 



(4.11) 



IC(s)|' = .5|CWI 




demonstrates how equation 4.4 has been extended to cover the range < 3?(s) < 1. 



It is easy to check the numerical validity of equation |4.11| by using the relationship 



shown in equation |2.1 
(4.12) 



IC(*)I'=IC(«)I»(^ 

Evaluating |C(s)|' at s = .1 + 2 0(K us ing equation 
evaluating |C( S )| using equation 



4.11 



4.12 



gives -28.52551836, while 
gives -28.52551645, where the sum over ze- 
ros (-0.55098684) is evaluated over the first 100000 zeros. A second example uses 
t = 9291.071149949, such that ((.5 + it) is close to a zero. For s = .1 + it, the sum 
over zeros is -3.8557529, with the two estimates being -199.4723 and -199.4719 re- 
spectively. Flipping to the other side, with s = .9+it, the sum over zeros is reversed 
(3.8557529), and the two e stimates are 0.2958258 and 0.2958254 respectively. 

and recall that 5ft ( — hi: ) = 2 ° 2 . 

y a-\-ib J a z +b z 



Consider equation 



4.11 
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(1) On the critical line, e quat ions 4.4 and 4.11 are identical, since at 5ft(s) = .5 



|x(s)| = 1 by Lemma 



and 



jft = Q since SJ(s) = ffi(p) = .5. 



(2) The sign of |C( S )| m equation 4.11 in the range < 5ft(s) < .5 depends 



on the sign of the terms in the brackets. The first term, 



|X(«)I' 



3.5 



is always 
The sign 



negative, since |x( s )l' is negative and monotonic by Lemma 
of the second term, involving the sum over the £ zeros, depends on the 
validity of the RH. For a single zero, with = .5, the contribution is 

negative since in this range, 5R(s — p) is negative. For a symmetric pair of 
zeros (supposing RH were false), say ph and pa, the contribution will be 
negative only if s is to the left of p^. If s is to the right of pr, and to the left 
of the critical line, that is 3?(pl) < 3t(s) < .5, then the total contribution 
of the pair is positive. 

Thus the RH cannot be established from a consideration of just the 
Hadamard product. □ 



Conjecture 4.11. 

process. 



The zeros of the £ function are derived from a single stochastic 



Since lim 



»(«)- 



, C(s) = 1, it follows directly from the functional equation 1.3 that 
: x( s )- Equation 4.11 demonstrates that, with the exception of 



limsR( s )^-oo C(s) 

the region around a £ zero, there is a high correlation between the behaviour of the 
£ function and the x function. The critical difference between the two is that by 
Lemma [3. 4| x( s ) has no zeros in S. We know that there are no £ zeros to the left 
of 5R(s) = 0, so the question is whether the ( function has become sufficiently close 
to the x function to the left of 3t(s) = .5 such that there are no zeros to the left of 
the critical line. 

We can get an indication as to whether the RH is true from a consideration of 



Lemma 4.7 We know from the work on random matrix theory that the distribution 



of C, zeros on the imaginary axis is a statistical phenomena. Thus, should these zeros 
also occur as symmetric pairs, then this too will also be a statistical phenomena. 
Let the nth £ zero, p n — a n + it n , and consider the statistic: 



(4.13) 
(4.14) 



p 



3? 



ti 



if 



Hence, if the RH is false, then J2 P ^ (p) consists of terms that are the ratio of two 
statistically derived variates - a n and t n . If we consider some standard statistical 
distributions that are similarly functions of the ratio of stochastic processes - beta 
(ratio of two gamma processes), F (ratio of two y 2 processes), Student's t (ratio of 
normal and y/ix 2 ) processes) - we observe that the distribution function cannot be 
expressed in an an alytic form, although the converse is not necessarily true. 

1.7 we note that J2 P ^ (p) = — log(2) — .51og(7r) + l + .57 where p are 
zeros of the C function, and 7 is Eulcr's constant. The fact that this 
sum can be expressed as an analytic form solely in terms of mathematical constants 
- and this is of course the crux of the conjecture - would imply that 3? fffj 



By Lemma 
the non trivia 
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does not consist of two stochastic processes. That would of course be true only if 
<r n = 0.5 V n - which would give us the RH. 
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